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ABSTRACT 

Context. PAMELA and, more recently, AMS-02, are ushering us into a new era of greatly reduced statistical uncertainties in exper¬ 
imental measurements of cosmic-ray fluxes. In particular, new determinations of traditional diagnostic tools such as the boron-to- 
carbon ratio (B/C) are expected to signiflcantly reduce errors on cosmic-ray diffusion parameters, with important implications for 
astroparticle physics, ranging from inferring primary source spectra to indirect dark matter searches. 

Aims. It is timely to stress, however, that the conclusions obtained crucially depend on the framework in which the data are interpreted 
as well as on some nuclear input parameters. We aim at assessing the theoretical uncertainties affecting the outcome, with models as 
simple as possible while still retaining the key dependencies. 

Methods. We compared different semi-analytical, two-zone model descriptions of cosmic-ray transport in the Galaxy: infinite 
slab(lD), cylindrical symmetry(2D) with homogeneous sources, cylindrical symmetry(2D) with inhomogeneous source distribu¬ 
tion. We tested for the effect of a primary source contamination in the boron flux by parametrically altering its flux. We also tested 
for nuclear cross-section uncertainties. All hypotheses were compared by minimisation techniques to preliminary results from 
AMS-02. 

Results. We And that the main theoretical bias on the determination of the diffusion coefficient index 6 (up to a factor two) is repre¬ 
sented by the assumption that no injection of boron takes place at the source. The next most important uncertainty is represented by 
cross-section uncertainties, which reach +20% in 6. As a comparison, nuclear uncertainties are more important than the shift in the 
best-fit when introducing a convective wind of velocity <30 km/s, with respect to a pure diffusive baseline model. Perhaps surpris¬ 
ingly, homogeneous ID vs. 2D performances are similar in determining diffusion parameters. An inhomogeneous source distribution 
marginally alters the central value of the diffusion coefficient normalisation (at the 10%, 1 cr level). However, the index of the diffusion 
coefficient 6 is basically unaltered, as well as the goodness of fit. 

Conclusions. Our study suggests that, differently for instance from the leptonic case, realistic modelling of the geometry of the Galaxy 
and of the source distribution are of minor importance to correctly reproduce B/C data at high energies and thus, to a large extent, 
for the extraction of diffusion parameters. The ansatz on the lack of primary injection of boron represents the most serious bias and 
requires multi-messenger studies to be addressed. If this uncertainty could be lifted, nuclear uncertainties would still represent a seri¬ 
ous concern; they degrade the systematic error on the inferred parameters to the 20% level, or three times the estimated experimental 
sensitivity. To reduce this, a new nuclear cross-section measurement campaign might be required. 
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1. Introduction 


The pattern of relative abundances of nuclei in the cosmic ra¬ 
diation (CR) is roughly similar to the one of the solar system 
material, with some notable exceptions: fragile nuclei (with low 
binding energies) such as or Li-Be-B are over-represented 
in CR. This CR component is usually interpreted as the result 
of production by spallation of heavier species during the prop¬ 
agation of primary cosmic rays—whose injected abundance is 
assumed to closely trace that of the solar system —in the inter¬ 
stellar medium. The ratios of these secondary to primary fluxes 
have long been recognised as a tool for constraining CR propaga¬ 
tion parameters, for some review see for instance |Maurin et al.| 
( 2002| ) and [Strong et al.|pOQ7| ). The boron-to-carbon ratio, or 
B/C, represents the most notable example among them. The key 
constraints on the diffusion parameters are inferred by its mea¬ 
surement, with the corresponding confidence levels (see for in¬ 
stance [MauHneT^ ([20^) widely used as benchmarks. It has 
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also been recognised that datasets available one decade ago were 
still insufficient for a sat isfactory measurement of these parame¬ 
ters ( [Maurin et al.|2QQl| ). 

The current decade is undergoing a major shift, however, 
with experiments such as PAMELA (Adriani et al. 2014) and 
most notably AMS-02 (http://www.ams02.org), which are 
characterised by significantly increased precision and better con¬ 
trol of systematics. The current and forthcoming availability 
of high-quality data prompts the question of how best to ex¬ 
ploit them to extract meaningful (astro)physical information. 
This new situation demands reassessing theoretical uncertain¬ 
ties, which will probably be the limiting factor in the param¬ 
eter extraction accuracy. As a preliminary work, preceding the 
actual data analysis, we revisit this issue to determine the rela¬ 
tive importance of various effects: some have already been con¬ 
sidered in the past, some were apparently never quantified. We 
also found that the main theoretical biases or errors are related 
to phenomena that can be described in a very simple ID diffu¬ 
sive model. We thus adopt it as a benchmark for our description. 
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reporting the key formulae that thus have a pedagogical useful¬ 
ness, too. In fact, we focus on determining the diffusion coeffi¬ 
cient, which we parameterise as conventionally in the literature 
(see for example Ptuskin et al. (1997) ): 


D(n) = Dofi 




^0 = 1 GV 


( 1 ) 


where Dq and S are determined by the level and power-spectrum 
of hydromagnetic turbulences, R is the rigidity, and the veloc¬ 
ity yS = v/c ^ 1 in the high-energy regime of interest here (ki¬ 
netic energy/nucleon > lOGeV/nuc). In fact, at lower energies 
numerous effects, in principle of similar magnitude, are present, 
such as convective winds, reacceleration, and collisional losses. 
At high energy, there is a common consensus that only diffusion 
and source-related effects are important. We focus on the high- 
energy region since it is the cleanest to extract diffusion param¬ 
eters, that is the least subject to parameter degeneracies. While 
adding lower-energy data can lead to better constraints from a 
statistical point of view, the model dependence cannot but grow. 
Since our purpose is to compare theoretical with statistical un¬ 
certainties from observations, our choice is thus conservative: in 
a global analysis, the weight of the former with respect to the 
latter is probably larger. To deal with a realistic level of statisti¬ 
cal errors of the data that will be available for the forthcoming 
analyses, we base our analyses on preliminary AM S-02 data of 
the B/C ratio ( |33rd Intern. Cosmic Ray Conf. 2Q13| ). 

This paper is organised as follows. In Sect.]2 we recall a sim¬ 
ple ID diffusion model providing our benchmark for the follow¬ 
ing analyses. This model certainly has pedagogical value, since it 
allows encoding the main dependences of the B/C ratio on input 
as well as astrophysical parameters in simple analytical formulae 
. At the same time, it provides a realistic description of the data, 
at least if one limits the analysis to sufficiently high energies. 
Relevant formulae are introduced in Sect. |2.1| while in subsec¬ 
tion [2^ we recall the main statistical tools used for the analysis. 
In Sect. we describe the main degeneracy affecting the analy¬ 
sis: the one with possible injection of boron nuclei at the sources. 
The next most important source of error is associated to cross- 
section uncertainties, to which we devote Sect. g In Sect.0 we 
discuss relatively minor effects linked to modelling of the geom¬ 
etry of the diffusion volume, source distribution, or the presence 
of convective winds. In Sect. we report our conclusions. 


2. B/C fit with a 1D modei 

2.1. ID diffusion model 


2H 


Free escape boundary 


Diffusive halo 


Infinitely thin matter disk 

Diffusive halo 


1 2h 


Free escape boundary 


Fig. 1. Sketch of the ID slab model of the Galaxy, with matter homoge¬ 
neously distributed inside an infinite plane of thickness 2h sandwiched 
between two thick diffusive layers of thickness 2H. 
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The simplest approach to model the transport of cosmic-ray 
nuclei inside the Calaxy is to assume that their production is 
confined inside an infinite plane of thickness 2h, that is sand¬ 
wiched inside an infinite diffusion volume of thickness 2H, sym¬ 
metric above and below the plane. The former region stands for 
the Galactic disk, which comprises the gas and the massive stars 
of the Milky Way, whereas the latter domain represents its mag¬ 
netic halo. A sketch of this model is given in Fig.[2 The bound¬ 
ary conditions fix the density of cosmic rays at the halo edges 
z = ±H to zero, while the condition h H (in practice, h is 
almost two orders of magnitude smaller than H) allows us to 
model the Galactic matter distribution as an infinitely thin disk 
whose vertical distribution is accounted for by the Dirac func¬ 
tion 2hS(z). Our focus on energies above lOGeV/nuc allows us 
to neglect continuous (ionisation and Coulomb) energy losses, 
electronic capture, and reacceleration. These subleading effects 
cannot be truly considered as theoretical uncertainties, since they 
can be introduced by a suitable upgrade of the model. However, 
taking them into account at this stage would imply a significant 
loss in simplicity and transparency. 

The well-known propagation equation for the (isotropic part 
of the gyrophase-averaged) phase space density if/a of a stable 
nucleus a, with charge (atomic number) Z^, expressed in units of 
particles cm“^ (GeV/nuc)“\ takes the form 


dt 




D 


dll/g 

dz 


Zmax 

2hS(z) • qa + d(z) ^ (Tb^a 

Zb>Za 

p 

-Siz) ■ (Ta ■ V - ij/a, 

^ISM 


p 

V - l]/b 

^ISM 


( 2 ) 


where the spatial diffusion coefficient D has been defined in 
Eq. 0. The cross-section for the production of the species a 
from the species b through its interactions with the interstel¬ 
lar medium (ISM) is denoted by whereas is the total 

inelastic interaction (destruction) cross-section of the species a 
with the ISM. The fragmentation of the nucleus b takes place at 
constant energy per nucleon, v hence stands for the velocities of 
both parent (b) and child (a) nuclei. The surface density of the 
Galactic disk is denoted by p, while misu is the average mass 
of the atomic gas that it contains. The values of the production 
cross-sections crb^a were calculated with the most recent formu¬ 
lae from Webber et al. (2003) The destruction cross-sections cr^ 
were computed by the semi-empirical formulae of Tripathi et al. 
(1997, 1999) The high-energy shapes of both cross-sections ex¬ 
hibit a plateau that allows one to approximate them as constants 
in this energy range. 

Solving the propagation Eq. 0 in the steady-state regime 
allows expressing the flux = (v/47r) if/a of the stable nucleus a 
inside the Galactic disk (z = 0) as 


Ja(Ek) = 


( ^max 

( Zb>Za 



where 


2D misM 
pvH 



+ 0-a}, 


(3) 


The fluxes J't of the parent species are also taken at z = 0. Here 
Qa, which stands for the source term, is homogeneous to a flux 
times a surface and is expressed in units of particles (GeV/nuc)“^ 
s“^ sr“^. It is related to qa through 


Qa 


i_^ =N i ^ 

in msM “llGV 


(4) 
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where Na is a normalisation constant that depends on the isotope 
a. We assumed an injection spectrum with the same spectral in¬ 
dex a for all nuclei. The value of Na should be adjusted by fit¬ 
ting the corresponding flux J'a to the measurements performed 
at Earth. However, these scarcely contain information on the iso¬ 
topic composition of cosmic rays. Nuclei with the same charge Z 
are in general collected together, irrespective of their mass. More 
isotopic observations would be necessary to set the values of the 
coefficients Na for the various isotopes a of the same element. 

In our analysis, we assumed solar system values (Lodders 
2003) for the isotopic fractions fa of the stable species a that 
were injected at the sources. We then proceeded by comput¬ 
ing the flux ffz of each element Z at Earth. We fixed the nor¬ 
malisation Nz for the total injection of all stable isotopes of 
the same charge Z by fitting the measured flux of that element. 
The normalisation entering in the calculation of Qa is given by 

= fa ' the sum of the fractions fa corresponding to the 
same element Z amounting to 1. The actual isotopic composition 
of the material accelerated at the source might be different from 
the solar system one, as is the case for neor (Binns et al. 2008) 
Our method might introduce a theoretical bias in CR element 
flux calculations. However, our main focus here is to extract the 
propagation parameters thanks to the different sensitivities be¬ 
tween primary carbon and (a priori) secondary boron. The only 
isotopes that come into play in the B/C ratio are the stable nu¬ 
clei ^^B, and ^^B, unstable plays a very minor role. 

The isotopes of either carbon or boron have similar rigidities and 
destruction cross-sections. Varying the isotopic composition of 
carbon (and of boron, should it be partially primary) does not af¬ 
fect the ratio calculation. Euthermore, secondary boron is mainly 
produced by the fragmentation of one particular isotope of each 
heavier element. Eor example, the primary component of is 
two orders of magnitude larger than that of This reduces 
the differences arising from the boron production cross-sections. 
Althought most of the isotopes at stake are stable, radioactive 
nuclei were also taken into account in the calculation, and we 
obtained more complicated expressions for the fluxes, which are 
not displayed here for brevity. They are reported for instance in 
Appendix A of |Putze et al.| ( |2010| ). By defining the total flux of 
a nucleus of charge Z as the sum over all its isotopes a 

Jz= Ja, ( 5 ) 

isotopes a 

of same Z 

and considering only the dominant contribution from stable nu¬ 
clei, the B/C flux ratio can be written as 


Jb(£^) 

JciEk) 



Jc 


/ jo-"-I- o-b) . (6) 


If we assume that there are no primary boron sources, that is, 
2 b = 0, this expression simplifies into 


fTBjEk) _ CTc^B CTb^B Jb .rj. 

Jc(Ek) + O-B 0-^'^ + 0 -B Jc 

The impact of relaxing this hypothesis is explored in Sect. 
where the effect of a non-vanishing value for 2 b is considered. 


2.2. Fitting procedure and benchmark values for this study 

We used the AMS-02 recent release of the B/C ratio ( |3 3rd Intern. | 
Cosmic Ray Conf.|2013| ) to study the impact of systematics on 


the propagation parameters. As explained above, we limited our¬ 
selves to the high-energy sub-sample, above 10 GeV/nuc. The 
set of Eqs. 0 is of triangular form. The heaviest element con¬ 
sidered in the network, which in our case is ^^Ee, can only suffer 
destruction. No other heavier species h enters in the determina¬ 
tion of its flux Ja, which hence is proportional to the injection 
term Qa. Once solved for it, the algebraic relation yields the so¬ 
lution for the lighter nuclei, down to boron. We evaluated the 
cascade down to beryllium to take into account its radioactive 
decay into boron. 

The primary purpose of our analysis is to determine the dif¬ 
fusion parameters Dq and d from the B/C flux ratio T = Jb /Jc- 
Another parameter of the model is the magnetic halo thickness 
H. As shown in Eq. Dq and H are completely degenerate 
when only considering stable nuclei, which provide the bulk of 
cosmic rays. In the following, H is therefore fixed at 4kpc for 
simplicity, although it should be kept in mind that, to a large ex¬ 
tent, variations in Dq can be traded for variations in H. Einally, 
the injection spectral index a also enters in the calculation of the 
B/C ratio through the source terms 2a- How strong its effect is 
on the best-fit diffusion parameters Dq and S is one of the ques¬ 
tions we treat in this section. To this purpose, we carried out a 
chi-square (x^) analysis of the B/C observations and minimised 
the function 



where the sum runs over the data points i whose kinetic ener¬ 
gies per nucleon are Ekj, while and ct/ stand for the cen¬ 
tral values and errors of the measurements. The theoretical ex¬ 
pectations also depend on the normalisation constants Na, 
which come into play in the source terms 2a of the cascade rela¬ 
tions 0- To determine them, we first fixed the spectral index a 
and the diffusion parameters Dq and 6. We then carried out an in¬ 
dependent;^^-based fit on the fluxes Jz of the various elements 
that belong to the chain that reaches from iron to beryllium. The 
measured fluxes are borrowed from the cosmic-ray database oi* 
Maurin et al. (2014) from which we selected the points above 10 
GeV/nuc. As explained above, this method yields the constants 
Nz and eventually the values of Na once the solar system iso¬ 
topic fractions fa are taken into account. The overall procedure 
amounts to profile over the normalisation constants Na to derive 
x\iQ as a function of a, 6 and Dq. Minimisations were performed 
by MINUIT (http: //www. cern. ch/minuit), a package inter¬ 
faced in the ROOT programme (https: //root. cern. ch). 

To check the accuracy and robustness of our fitting proce¬ 
dure, a preliminary test is in order. A commonly accepted no¬ 
tion is that the B/C ratio does not depend, to leading order, on 
the spectral index a. There is indeed no dependence on a in the 
cross-section ratios of Eq. (|7]) in the pure diffusive regime where 
(Tb . We have checked numerically that this behaviour 

holds by calculating the B/C best-fit values of the diffusion pa¬ 
rameters at fixed spectral index a. The results are reported in 
Eig.[^ where Dq and 6 are plotted, with their confidence limits, 
as a function of a. We scanned over the physical range that ex¬ 
tends from -2.5 to -2 and observed that the relative variations 
of Dq and ^ are 5% and 1%, respectively. The blue {Dq) and red 
(b) bands are almost horizontal. An anti-correlation between Dq 
and 6 is marginally noticeable and can be understood by the in¬ 
terplay of these parameters inside the diffusion coefficient D, the 
only relevant parameter that the B/C fit probes. We attribute the 
small variation of Dq with a to the different sensitivities of the 
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Fig. 2. Relative variations of the best-fit propagation parameters (com¬ 
pared to the benchmark model of Tablewith respect to the injection 
spectral index a. 


Fig. 3. Preliminary AMS-02 measurements of the B/C ratio ( |33rd In- 
|tern. Cosmic Ray Conf.|2013] ) are plotted as a function of kinetic energy 
per nucleon. The theoretical prediction of the ID/slab reference model 
of Tab.[2is also featured for comparison. 


normalisation constants Nz of nitrogen and oxygen to the low- 
energy data points as compared to carbon. This could result in 
fluctuations of the N^/Nq and Nq/Nq ratios with respect to the 
actual values. In any case, the extremely small dependence of the 
B/C ratio on a confirms the naive expectations and suggests that 
it is useless and simply impractical to keep cr as a free parameter. 

Nonetheless, there is a particular value of the injection index 
that best fits the fluxes of the elements Z that come into play 
in the cascade from iron to beryllium. By minimising the 
function 



^z,i ] 


(9) 


we find a = -2.34 as our benchmark value. Applying then our 
B/C analysis yields the propagation parameters Dq and d of the 
reference model of Table[2 which we used for the following anal¬ 
yses. The corresponding B/C ratio is plotted in Fig. as a func¬ 
tion of kinetic energy per nucleon and compared to the prelim¬ 
inary AMS-02 measurements ( |33rd Intern. Cosmic Ray Conf.| 
|2013| ). In what follows, we study how and 6 are affected by 
a few effects under scrutiny and gauge the magnitude of their 
changes with respect to the reference model. We could have de¬ 
cided to keep the injection index a equal to its fiducial value of 
-2.34, but we preferred to fix the spectral index y = a-d = -2.78 
of the high-energy fluxes J'z at Earth. Keeping a fixed would 
have little effect on the B/C ratio, but would degrade the good¬ 
ness of the fits on absolute fluxes. 


Reference parameter values 

a 

-2.34 

Do [kpe^/Myr] 

(5.8 + 0.7) • 10-2 

d 

0.44 ± 0.03 


5.4/8 0.68 

y = a - S (fixed) 

-2.78 


Another crucial test of our fitting procedure is to check how 
the results depend on the low-energy cut-off E^xxi above which 
we carried out our analysis. We set the flux spectral index y to 
its benchmark value of Table [T] and determined the B/C best-fit 
values of the diffusion parameters as a function of which 
was varied from 5 to 30GeV/nuc. The results are plotted in 
Fig. 1^ with the Icr and 2(t uncertainty bands. As expected, 
the statistical errors increase when moving from a low £’cut to 
a higher value. That is why the reduced (dashed line) de¬ 
creases steadily as the cut-off energy is increased. The higher the 
cosmic-ray energy, the fainter the flux and the scarcer the events 
in the detector. The widths of the blue (Dq) and red (^) bands 
at E’en! = lOGeV/nuc, however, are not significantly larger than 
for a cut-off energy of 5 GeV/nuc. This suggests that our esti¬ 
mates for the statistical errors are slightly pessimistic, which is 
acceptable and consistent with our purpose. 

The other trend that we observe in Fig. is a shift in the 
preferred value of 6 to increasingly lower values as we limit the 
analysis to increasingly higher energies. This is no limitation of 
our procedure. On the contrary, it is a real feature that the data 
exhibit, as is clear in Fig.|^ where the tail of the B/C points does 
look flatter above 50 GeV/nuc. The anti-correlation between 6 
and Dq that we observe in Fig. has already been explained by 
the interplay of these two parameters inside the diffusion coef¬ 
ficient D, to which the B/C ratio is sensitive. The increase of 
Dq is then generic and does not signal any new effect. At that 
stage, the statistical uncertainties are still of the same order as 
the systematic uncertainties generated by using different energy 
cuts. Should the decrease of S with Ecut be confirmed with higher 
statistics, some intrinsic explanation might be necessary for the 
failure of a power-law fit. See for instance Sect. for a possible 
explanation. 

3. Primary boron? 


Table 1. Benchmark best-fit parameters of the ID/slab model, with re¬ 
spect to which comparisons are subsequently made. 


Typical fits of the B/C ratio are based on the assumption that 
no boron is accelerated at the source, so that the term propor¬ 
tional to 2 b at the right-hand side of Eq. ([^ vanishes. How¬ 
ever, this is just an assumption that need to be tested empiri¬ 
cally. It is crucially linked to the hypothesis that the accelera- 
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Fig. 4. Relative variations of the best-fit propagation parameters, as 
compared to the benchmark model of Tab. with respect to the low- 
energy cut-off Ecxxi above which we carry out the B/C analysis. 


tion time is much shoter than the propagation time within the 
magnetic halo and that it occurs in a low-density environment. 
On the other hand, typical astrophysical accelerators such as su¬ 
pernova remnants might have the capability to accelerate up to 
TeV energies for of about ~ 10^ years in an interstellar 
medium with ^ism ~ 1 cm“^, or greater when surrounded by 
denser circumstellar material. The corresponding surface density 



z 

5 

TY. 


[a 

Do 


lij 





















2 4 6 8 10 12 

Primary abundances ratio N^q/Nq [%] 


^iSM c fiife ~ 10^^ cm“^ easily leads to percent-level probabilities 
for nuclei to undergo spallation in the sources. A factor of only a 
few times higher than this would certainly have dramatic conse¬ 
quences on the information inferred from secondary-to-primary 
ratios. More elaborate versions of this idea and related phe¬ 
nomenology have also been detailed as a possible explanation of 
the hard spectrum of secondary positron data ( |Blasi|20()9} |Blasi| 
|& Serpic^|2009t [Mertsch & Sar kar||2009|), which was recently 
compared with the AMS-02 dats (Mertsch & Sarkar 2014) 

Apparently little attention has been paid to the bias intro¬ 
duced by the ansatz 2 b = 0- To the best of our knowledge, we 
quantify it here for the first time. As can be inferred from Eq. 
in the presence of a primary source 2 b, the B/C ratio exhibits 
a plateau as soon as the cross-section ratio + ctr) 

becomes negligible with respect to the primary abundances ra¬ 
tio Ab/A^c- This happens at sufficiently high energy since 
increases with the diffusion coefficient D. The height of this 
high-energy B/C plateau is approximately given by the value of 
Ab/Ac. In the presence of this behaviour, the spectral index 6 
must increase to keep fitting the data at low energy, that is, here 
around lOGeV/nuc. This also implies that Dq decreases with 
Ab / Ac as a result of the above-mentioned anti-correlation be¬ 
tween the diffusion parameters. 

We have thus scanned the boron-to-carbon ratio at the source 
and studied the variations of the best-fit values of and d with 
respect to the reference model of Table Our results are illus¬ 
trated in Fig.|^ where the left panel features the confidence lev¬ 
els for d (red) and Do (blue) as a function of the Ab /Ac ratio. The 
B/C fit is particularly sensitive to the last few AMS-02 points, 
notably the penultimate data point, around 214 GeV/nuc, for 
which the B/C ratio is found to be ~ 9%. In the right panel, the 
theoretical expectation for that point is plotted (solid red curve) 
as a function of the primary abundances ratio, while the dashed 
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Fig. 5. Left panel: variations of the best-fit propagation parameters 
Dq (blue) and 5 (red) relative to the benchmark values of Table [T] as 
a function of the primary boron-to-carbon injection ratio. The refer¬ 
ence model corresponds to the conventional no boron hypothesis for 
which Ab/Ac vanishes. Right panel: the theoretical value of the B/C 
ratio at 214 GeV/nuc (solid red curve) is plotted as a function of the pri¬ 
mary boron-to-carbon injection ratio. The dashed black curve indicates 
the goodness of the B/C fit. As long as Ab/Ac does not exceed 13%, 
the theoretical B/C ratio is within 2cr from the AMS-02 measurement 
(dashed-dotted green curve). 


black curve indicates how the goodness of fit varies. It is interest¬ 
ing to note that a minor preference is shown for a non-vanishing 
fraction of primary boron, around 8%, due to the marginal pref¬ 
erence for a flattening of the ratio already mentioned in the pre¬ 
vious section. The Ab / Ac ratio is only loosely constrained to be 
below 13%. Such a loose constraint would nominally mean that a 
spectral index d more than three times larger than its benchmark 
value would be allowed, with a coefficient Do one order of mag¬ 
nitude smaller than indicated in Table In fact, such changes 
are so extreme that they would clash with other phenomenologi¬ 
cal or theoretical constraints and should probably be considered 
as unphysical. A spectral index 6 in excess of 0.9, correspond¬ 
ing to a relative increase of 100% with respect to our benchmark 
model, is already so difficult to reconcile with the power-law 
spectrum of nuclei and the present acceleration schemes that it 
would probably be excluded. The message is quite remarkable 
however. The degeneracy of the diffusion parameters with a pos¬ 
sible admixture of primary boron is so strong that it dramatically 
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degrades our capability of determining the best-fit values of Dq 
and S, and beyond them the properties of turbulence, unless other 
priors are imposed. 


4. Cross-section modelling 

The outcome of cosmic-ray propagation strongly depends on the 
values of the nuclear production and destruction cross- 
sections with the ISM species, mainly protons and helium nuclei. 
Some of these are measured, albeit in a limited dynamical range, 
while a significant number of them rely on relatively old semi- 
empirical formulas, calibrated to the few available data points. In 
this section, we discuss how parametric changes in these inputs 
reflect on the B/C ratio. The effect of cross-section systematics 
was already studied by Maurin et al. (2010), who parameterised 
it in terms of a systematic shift with respect to the energy. Since 
we consider here only the high-energy limit, we simply allowed 
for a rescaling of the cross-sections. However, we distinguished 
between two cases: a correlated (/'/') or anti-correlated (/'\) 
rescaling between the production and the destruction 
cross-sections. These in fact are not affected by the same uncer¬ 
tainties. It is often the case that the latter are known to a bet¬ 
ter precision then the former since they rely on a richer set of 
data. A priori, it is conceivable that several relevant production 
cross-sections might be varied independently. It is worth noting, 
however, that only a few nuclei - notably oxygen and carbon 
(~ 80%), and to a lesser extent nitrogen (~ 7%) - are in fact 
responsible for most of the produced boron, as shown in Fig.|^ 
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nik et al. 1998) (see Moskalenko et al. (2001) and Moskalenko 
& Mashnik (2003) for a more complete list of references). In this 
work we compare the values given directly by the default cross- 
section parameterisations without any renormalisation (which 
can be implemented however). 

In the case of the destruction cross-sections cr^, we compared 
our reference model (Tripathi et al. 1997) with the parameteri¬ 
sations of |Barashenkov & Polanski] (19941 ), Letaw et al. (1983) 
and Wellisch & Axen (1996) The last case only applies to el¬ 
ements with Z > 5, while the Letaw et al. (1983) modelling is 
conserved for lighter nuclei. Figure shows the relative differ¬ 
ences between our reference model and the three other semi- 
empirical approaches and allows deriving an indicative lower 
limit on the systematic uncertainties for the destruction cross- 
sections of roughly 2 to 10% for the B/C ratio. The systematic 
difference is at the 3% level for the channels (CNO) that con¬ 
tribute most to secondary boron production. The difference to 
our reference model is stronger for larger charges (Z > 10), but 
these nuclei have a negligible contribution to the B/C ratio. 



Fig. 7. Relative differences between our reference mode] (Tripathi et al. 
1997) for the destruction cross-sections cr^ and the other parameteri¬ 
sations by Letaw 1983 (Leta w et al. 1983) Wellish 1996 (We llisch & 
Axen 1996) and B&P 1994 ( [Barashenkov & Polanski|p^^ are dis¬ 
played as a function of the nucleus charge, at an energy of 10 GeV/nuc. 
Each bin is characterised by a given charge Z and encodes the arithmetic 
mean over the corresponding isotopes. Only the elements involved in 
the cascade from iron to beryllium are displayed. 


Fig. 6. Contribution of the various primary nuclear species to the sec¬ 
ondary boron flux at 10 GeV/nuc, as estimated with the semi-empirical 
code by Webber 03. 

First, we need to assess the reasonable range over which the 
various cross-sections of the problem are expected to vary. For 
this, we compared our reference models for the destruction and 
production cross-sections with those used in popular numerical 
propagation codes such as GALPROF gtrong & Moskalenko 
2001) and DRAGON (Evoli et al. 2008|] The database imple¬ 
mented these two codes traces back to the GALPROP team work 
and is based on a number of references including - but not lim¬ 
ited to - Nuclear Data Sheets and Los Alamos database (Mash- 

^ Updated version of these two codes can be found 
at: https: //sourceforge .net/pro jects/galprop and 

http: //WWW . dragonproj ect. org/Home. html respectively. 
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For the production cross-sections one may chose be¬ 

tween the semi-empirical approach proposed b> Silberberg et al. 
(1998), subsequently revised in 2000 and called here S&T 00, 
and the parameterisation provided b> Webber et al. (1990) (here¬ 
after Webber 93) and its updates of 1998 (Webber et al. 1998) 
and 2003 (Webber et al. 2003) We selected the last set of values 
as our reference model, to which we have compared the other pa¬ 
rameterisations to gauge the uncertainties that affect, on average, 
the values of (Tb^a- The relative differences between Webber 93, 
Webber 98, and S&T 00 with respect to Webber 03 are plotted 
in the form of the three histograms of Fig. The charges of the 
parent and child nuclei are given on the vertical and horizon¬ 
tal axes. The most important reactions, whose cross-sections are 
higher, correspond to a change of charge AZ not in excess of 3 
during the fragmentation process and are located close to the di¬ 
agonals of the 2D-grids of Fig. We first note that the Webber 
93 and 98 production cross-sections are on average larger than 
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Fig. 8. 2D histograms feature the relative differences between various 
semi-empirical models currently used to calculate the production cross- 
sections cTb^a’ Our reference model is Webber et al. (2003) (Webber 
03), and we compare it to the parameterisations fronc Webber et al. 
(1990) (Webber 93), Webber et al. (1998) (Webber 98) anc Silberberg 
et al. (1998) (S&T 00). The charges of the parent and child nuclei are 
given on the vertical and horizontal axes, respectively. The relative dif¬ 
ference in each bin is given by the arithmetic mean over the various 
isotopes of each element. A detailed view provides the most important 
channels for the B/C ratio studies. For a fragmentation of AZ < 4, we 
also give the first and second moments of the uncertainty distributions. 


the values of the Webber 03 reference model. Most of the pixels 
on the diagonals of the corresponding histograms are red, and we 
measured an excess yu on the reactions for which AZ < 4 of 18% 
and 9.7% for Webber 93 and Webber 98 as compared to Webber 
03. Furthermore, in both cases the dispersion of these differences 
is quite large and amounts to 31% for Webber 93 and 30% for 
Webber 98. A rapid comparison between S&T 00 and Webber 03 
would also leave the impression that in the former case, the reac¬ 
tions in the upper left comer of the histogram have cross-sections 
exceedingly larger than for the Webber 03 parameterisation. A 
close inspection along the diagonal indicates, on the contrary, 
that the S&T 00 values for AZ < 4 are on average 13% higher 
than for the reference model, with a dispersion cr of 28% similar 
to the other cases. The main production channels of secondary 
boron are listed in Tablej^and are also displayed in the expanded 
views of the small square regions that sit in the lower left corners 
of the histograms of Fig.[^ The most relevant reactions involve 
the stable isotopes of carbon, nitrogen, and oxygen fragmenting 
into and and are indicated in boldface in TableThe 
largest contributor to secondary boron is The three semi- 
empirical models with which we compared our Webber 03 refer¬ 
ence parameterisation tend to predict production cross-sections 
that are 15% for S&T 00 to 25% for Webber 93 larger. In con¬ 
trast, those models underpredict the spallation of by 10% in 
the case of Webber 93 and 98 to 18% for S&T 00. In the latter 
case, the production cross-section of ^^B from is 68% larger 
than for Webber 03. But nitrogen only contributes ~ 7% of the 
secondary boron, and this has no significant impact. To sum¬ 
marise this discussion, the production cross-sections cr^^a can 
be varied up or down by a factor of order 10-20% with respect 
to Webber 03. 

Varying the various production and destmction cross- 
sections has an effect on the calculation of the B/C ratio and thus 
affects the determination of the propagation parameters Dq and S. 
Before gauging this effect, we remark that secondary boron is es¬ 
sentially produced by CNO nuclei, as indicated in Fig. These 
are essentially primary species for which J't is approximately 
given by the ratio Qb/icr^^^ + ctz?) and is proportional to the in¬ 
jection normalisation Nb. Furthermore, the relevant destmction 
cross-sections (Tc, ctn and (Tq being approximately equal to each 
other, with an effective value ranging from 290 to 317 mb, we 
conclude that the flux ratios J'bl.Jc are given by the correspond¬ 
ing ratios Nb/Nc of the injection normalisation constants, with 
the consequence that relation (|7]) simplifies to 

JBiEk) ^ y CTh^B ^ 

Jc(Ek) + Ac' ^ ^ 

As mentioned at the beginning of this section, we first rescaled in 
our code all production (Tb^a and destruction cross-sections 
by the same amount k, which ranges from 0 to 2, to study how 
Do and S are affected by this change. The results are summarised 
in the left panel of Fig. The diffusion index S does not suffer 
any change, whereas the diffusion normalisation Dq increases 
linearly with the rescaling factor k. Multiplying both and 
(Tb by the same factor k in Eq. ( p^ amounts to dividing the dif¬ 
fusion cross section by k. The B/C ratio depends then on the 
ratio Ik, which scales as Dq/k. The theoretical prediction on 
the B/C ratio is not altered as long as that ratio is kept constant, 
hence the exact scaling of Dq with k displayed in the left panel 
of Fig.|^ The energy behaviour of the B/C ratio is not sensitive 
to the rescaling factor k, which has been absorbed by Dq, and the 
fit yields the same spectral index S irrespective of how much the 
cross-sections have been changed. Despite the relatively modest 
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Table 2. Comparison between different cross-section parameterisations for the main production channels of secondary boron. The reference model 
used in our calculations of the fluxes is adapted fron: Webber et al. (2003) (Webber 03) and is compared to previous releases b> Webber et al. 
(1990) (Webber 93) and Webber et al. (1998) (Webber 98) as well as to the work fron: Silberberg et al. (1998) (S&T 00). The dominant production 
channels, which involve the stable isotopes of carbon, nitrogen, and oxygen, are listed in boldface. 


alterations, the effect discussed here has two qualitatively inter¬ 
esting consequences. To commence, a systematic uncertainty on 
the central value of Dq at the 5 to 10% level seems unavoidable 
due to the current uncertainty level of about 10% on the nuclear 
cross-sections. Then, fully correlated changes in both production 
and destruction cross-sections can break the degeneracy between 
Dq and S. 

We now analyse the effects of an anti-correlated change of 
the production and destruction cross-sections. Surpris¬ 

ingly, this has never been considered before, as far as we know, 
although the potential effect of this rescaling clearly is very 
strong. Multiplying by a factor k while rescaling ctb by 

a complementary factor of (2 - k) leads to the B/C ratio 

JuiEQ y f ] Nb 

JciEk) + Nc 

Keeping the B/C ratio constant while increasing a: at a given en¬ 
ergy translates into keeping the ratio 

+ 2o-b _ CE ^ + 2o-b 

K K 

roughly constant, where C is a constant directly proportional to 
Dq. It can be immediately inferred that, when k increases, C and 
Dq have to increase and thus d has to decrease. This trend is 
confirmed in the right panel of Fig.From realistic assessments 
of the minimum systematic uncertainties of about 10% derived 
from the different cross-section models, we estimate a systematic 
uncertainty of 10% on 6 and of 40% on Dq. 


5. Systematics related to CR propagation modelling 

A significant effort has been made in recent years to provide in¬ 
creasingly sophisticated modelling of the CR diffusion environ¬ 
ment, source distribution, and alternative forms of CR transport. 
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In this section we discuss a perhaps surprising conclusion: these 
effects are less relevant for the prediction of B/C than the effects 
discussed previously (which are instead usually neglected)! The 
message is: although the efforts invested by the community in 
refining CR propagation modelling could have and have had im¬ 
portant implications for other observables, for the mere purpose 
of fitting B/C to infer diffusion propagation parameters they are 
to a large extent unnecessary complications, until one can signif¬ 
icantly reduce the biases previously discussed. 


5.1. Geometric effects 


The crude modelling of the diffusive halo as an infinite slab may 
appear too simplistic. In this section, we estimate the effects of 
a 2D cylindrical diffusion box, modelled as in Fig. ( p^ . Fur¬ 
thermore, we assess the effect of adding a radial dependence in 
the injection term, as opposed to the uniform hypothesis. These 
can be seen as upper limits to reasonable systematics due to 
simplified description of the spatial dependence of the diffu¬ 
sion medium or source term: given our limited knowledge on 
this subject, even the most detailed modelling of the propagation 
medium and source term, in fact, may not be fully realistic. 


The formalism in such a situation is well known and we do 
not repeat it here (it has been summarised for instance in |Putze| 
et al.| ( |201^ ). It suffices to say that to take advantage of the cylin¬ 


drical symmetry, Eq.Q can be projected on the basis of the zero 
order Bessel functions /^(r) = Jq j ensuring that the den¬ 


sity vanishes on the edge of the cylinder of radius Rgai = 20kpc. 
The flux of an isotope is then the sum over all its harmonic com¬ 
ponents 


Ja(Ek,Ro) = J]jo(^l^^X 


a(Ek). 


(13) 
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Fig. 9. Effect of rescaling nuclear cross-sections for boron production 
channels and destruction ones: the left panel assumes correlated, the 
right panel anti-correlated rescaling. 



Fig. 10. Cylindric model: the matter is homogeneously distributed in¬ 
side a thin disk of thickness 2h and radius Rgai = 20kpc. The solar 
system is at - 8 kpc from the Galactic centre. 


The results, reported in Tableallow us to draw a few con¬ 
clusions: 

- the presence of a new escape surface at Rg^i - 20 kpc is ba¬ 
sically irrelevant: the best-fit S and its error remain the same, 


with a statistically insignificant, 2% modification of the best- 
fit value of Do ; 

- perhaps more surprisingly, even the replacement of a uni¬ 
form source distribution with a commonly assumed donut 
distribution of the forn: (Yusifov & Kiiguk 2004) 


q(r) oc 


I r-hO.55 
\Rq -f 0.55 


1.64 


exp -4.01 


r-Rp 
Rq -f 0.55 


(14) 


has minor effects, a mere 1% modification in the best-fit de¬ 
termination of S, and a ~ 13% lowering of the best-fit value 
of Do, still statistically insignificant (roughly a 1 cr effect); 

- since the goodness of fit is similar, the B/C observable is es¬ 
sentially insensitive to these improvements. Unless they are 
justified by the goal of matching or predicting other observ¬ 
ables, the complication brought by the 2D modelling of the 
problem are unnecessary in achieving a good description of 
the data. 


5.2. Convective wind 


Although the high-energy CR propagation is mostly diffusive, 
the advection outside the Galactic plane (for instance due to stel¬ 
lar winds) has a non-negligible effect, which we now quantify. 
We adopted the simplest model of constant velocity wind, di¬ 
rected outside the galactic plane, with magnitude u. Taking this 
effect into account this effect, the ID, stationary propagation 
equation can be written as 


^ /r. ^ / 'l , ^ ^ ^ (Idu iEk + 2m) 

- & ("s'*") M (3 

^max 

+ S{z)craV-^tl/a = 2hS{z)qa + d{z) Y CTb^aV-^lJ/b, (15) 


Zb>Za 


^ISM 


The two new terms (second and third one on the LHS) account 
for the advection of the cosmic-ray density and the adiabatic 
losses, respectively. A characteristic time of these two processes 
can be estimated inside the thin disk of matter : 


^ advection 


and 


h 
u 

= 5 • 


0.1 kpc 
20km/s 

0.1 kpc j \ 


1 20 km/ s 


My, 


^ adiabatic 


1 

M 

15- 


(Vt/) ^3A- 

/ u 

h \ / 20km/s \ 

0.1 kpc / \ / 


.(d)M. 


(16) 


(17) 


This means that adiabatic losses can be safely neglected com¬ 
pared to the typical diffusion time of 


^diffusion ("^ > 10 GV) < Tditfusion(lOGV) 

h 


hH 


= 2 - 


0.1 kpc 


t—) 

\4kpc/ 


D(IOGV) 

5.8-10-2 •(2-10)^44 kpcVMy 
D 


My. 


(18) 


It is clear that our previous results provide a suitable first-order 
approximation at least at high energy, with the leading correction 
at energies near 10 GeV/n given especially by the advection. The 
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Geometry 

Plane - ID 

Cylindrical - 2D 
homogeneous source distribution 

Cylindrical - 2D 
realistic source distribution 

Do [kpc^/Myr] 

(5.8 + 0.7) • 10-^ 

(5.7 ± 0.7) • 10“^ 

(5.0 + 0.6) • 10-^ 


N/A 

-2% 

-13% 

5 

0.441 +0.031 

0.439 + 0.031 

0.445 + 0.032 


N/A 

0% 

+ 1% 

4/c/ndof 

5.4/8 ^ 0.68 

5.4/8 ^ 0.68 

5.5/8 ^ 0.69 


Table 3. Results on the propagation parameters fitted on the B/C for different geometries. 


adiabatic energy loss, instead, is several times smaller and can be 
safely ignored in the following. 

The solution of Eq. © neglecting adiabatic losses has the 
same form of Eq. ^ for the flux of stable species, modulo the 
change 


D^D' = 


Hu 


.exp(-iif)' 


(19) 


SO that the behaviour of the solution smoothly interpolates be¬ 
tween the convective timescale at low energy and the diffusive 
one at high energy: this can be simply checked by neglecting the 
exponential with respect to unity for a high value of its argument, 
or Taylor-expanding it to first order in the opposite limit. This 
formula also suggests that, if one fits the data by neglecting the 
convective wind, one biases its result towards a lower value of S, 
and a corresponding higher value of D, so to reproduce a flatter 
dependence with energy at low-energy as for the case described 
by Eq. ( p^ , as illustrated in Eig. [TT] Quantitatively, a variation 
of 15 km/s in u is roughly similar to a Icr shift in the benchmark 
parameters. Note, however, that the goodness of the fit worsens, 
or in other words, high-energy data are better described by a pure 
diffusive behaviour than by a convective-diffusive one. Overall, 
we conclude that these effects appear still somewhat less impor¬ 
tant in determining the diffusion parameters from high-energy 
data than the role of primary boron or even cross-section uncer¬ 
tainties. While convection, adiabatic losses, reacceleation, etc. 
are important to account for when extending the analysis down 
to very low energies (sub-GeV/nuc) or in global analyses, they 
do not currently constitute the main limitations to the determina¬ 
tion of Dq or S from high-energy data. 


6. Conclusion 

The high-precision measurements of cosmic-ray fluxes that have 
become available in recent years prompt the question of the the¬ 
oretical uncertainties inherent to the models used to interpret 
them. In this article, we have compared the effect of different 
theoretical biases with statistical uncertainties in the determina¬ 
tion of diffusion parameters from the boron-to-carbon flux ratio, 
or B/C. This is representative of a much broader class of observ¬ 
ables, involving ratios of secondary to primary species, which 
have been recognised as key tools for diagnostics in cosmic-ray 
astrophysics. We adopted a pedagogical approach, showing and 
interpreting the results whenever possible within simple analyt¬ 
ical models. We also used preliminary AMS-02 data and lim¬ 
ited the analysis to energies above lOGeV/nuc, which gives a 
pessimistic—hence conservative—estimate of the statistical un¬ 
certainties that will eventually be available. 

Our main results, summarised in Tableare the following: 

- The single most important effect that we quantified (to the 
best of our knowledge, for the first time) is the degeneracy 
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Wind velocity (km/s) 


Fig. 11. Variations of the best-fit propagation parameters with respect 
to the velocity of the convective wind. 


between diffusion parameters and a small injection of pri¬ 
mary boron at the source, finding at present even a statisti¬ 
cally insignificant preference for a small but finite value for 
a primary boron flux. This degeneracy cannot be removed 
by high-precision measurements of B/C, but probably re¬ 
quires multi-messenger tests and certainly demands further 
investigations, in particular if data should manifest a signifi¬ 
cant preference for a high-energy flattening of secondary-to¬ 
primary ratios. 

- The second most important theoretical uncertainty is asso¬ 
ciated to cross-sections. In particular, anti-correlated modifi¬ 
cations in the destruction and production cross-sections with 
respect to reference values may also have an effect on the de¬ 
termination of the diffusion index S, another effect discussed 
here for the first time. This should be kept in mind when 
comparing the outcome of data analyses relying on differ¬ 
ent databases for cross-sections. The good news is that this 
problem is not due to intrinsic limitations in the astrophys- 
ical modelling or the lack of astrophysical data, but to the 
scarce laboratory measurements available. Eor the case of 
boron, experiments of production cross-sections via spalla¬ 
tion of oxygen, carbon and, to a minor extent, nitrogen, are 
essentially what would be needed to set the predictions on 
much firmer grounds. 

- Other effects we tested for are typically less important and 
are similar to or smaller than statistical uncertainties: effects 
such as those of convective winds, certainly important in 
more complete analyses including low-energy data, appear 
unlikely to bring uncertainties large enough to compete with 
the above-mentioned uncertainties. We also showed how the 
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Wind 

1D/2D geometry 

Cross-sections 

Primary boron 

ADqIDq 

-40% 

-2 to-13% 

+60% 

0 to -90% 

ASIS 

-fl5% 

0 to -fl% 

±20% 

0 to+100% 


Table 4. Summary of the main systematics found in current analyses in determining the propagation parameters by fitting the B/C ratio. 


geometry of the diffusive box and the distribution of sources 
is virtually irrelevant, at least if only a B/C data analysis 
is concerned. More or less realistic radial distribution of 
sources, while it may marginally affect the determination of 
Dq, is still indistinguishable from the goodness-of-fit point of 
view. Another outcome of this exercise is that at least at the 
10% level, Do is degenerate with a choice of geometry and 
source distribution, in addition to the already well-known de¬ 
generacy with the diffusive halo height H. 

In conclusion, we found that the main uncertainties in infer¬ 
ring diffusion parameters from B/C (and we expect from other 
secondary-to-primary ratios, too) depend on theoretical priors on 
sources (linked to sites and mechanisms of acceleration!) and, to 
a lesser extent, to nuclear cross-sections. While exploring more 
complicated schemes and geometries for the diffusion may thus 
be important, we can anticipate that sensitivity to such effects 
will probably require fixing more mundane questions first! A 
multi-messenger strategy, coupled to a new measurement cam¬ 
paign of nuclear cross-sections, appears to be a next crucial step 
in that direction. 
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